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0\ ■ Abstract 



^3 



We show that in some regions of supersymmetric parameter space, CP vi- 
olating effects that mix the CP-even and CP-odd Higgs bosons can enhance 



I. INTRODUCTION 



<^ ' the neutrahno annihilation rate, and hence the indirect detection rate of neu- 

tralino dark matter, by factors of 10^. The same CP violating effects can 
reduce the neutralino scattering rate off nucleons, and hence the direct detec- 
tion rate of neutralino dark matter, by factors of 10""^. We study the depen- 
^ ' dence of these effects on the phase of the trilinear coupling A, and find cases 

' in the region being probed by dark matter searches which are experimentally 

. excluded when CP is conserved but are allowed when CP is violated. 

00 
O 

a^ 

p • The nature of the dark matter in the universe is one of the outstanding questions in 

p • astro/particle physics. One of the favored candidates is the lightest supersymmetric (SUSY) 

D ! particle. Such a particle is weakly interacting and massive (with mass in the range 1 GeV - 

^ ' few TeV), and hence is frequently characterized as a WIMP (weakly interacting massive 
particle). In the minimal supersymmetric standard model (MSSM), the lightest SUSY 

^ ■ particle in most cases is the lightest neutralino, a linear combination of the supersymmetric 

■ partners of the photon, Z° boson, and neutral-Higgs bosons, 

X? = NuB + NuW' + NisHi + N^Jl^ (1) 

where B and W'^ are the supersymmetric partners of the U(l) gauge field B and of the third 
component of the SU(2) gauge field that mix to make the photon and Z*^ boson. (We 
will also use the letter x for Xv) 

Much work has been done studying the possibilities for detecting these particles. Possi- 
bilities include direct detection [Q, whereby the particle interacts with a nucleon in a low 
temperature detector, and is identified by the keV of energy it deposits in the nucleon; and 
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indirect detection, whereby (1) the particles are captured in the Sun or Earth, sink to the 
center of the Sun or Earth, annihilate with one another in the core, and give rise to particles 
including neutrinos which can be detected by experiments on the surface of the Earth 0, or 
(2) the particles annihilate in the galactic halo and produced anomalous components in the 
flux of cosmic rays . The interaction processes of the lightest SUSY particle are clearly of 
great importance in calculations of predicted rates for both direct and indirect detection. 

In this paper we discuss the effect of CP violation on the neutralino annihilation and 
scattering cross sections. The MSSM introduces several new phases in the theory which 
are absent in the standard model. Supplemented by a universality condition at the grand 
unification scale, only two of these are independent. In this case, one may choose to work 
in a basis in which the two non-trivial CP-violating phases reside in /i and the universal 
soft trilinear coupling A of the Higgs fields to the scalar fermions /. Previously Falk, Ferstl 
and Olive have considered the effect on neutralino cross sections of a nonzero phase of /i, 
the mixing mass parameter involving the two Higgs chiral superfields in the superpotential. 
Here, on the other hand, we consider the effect on neutralino cross sections of the case where 
the soft trilinear scalar couplings Af are all complex numbers, where subscript / refers to 
the quarks. To be specific, we will take A = At = Ab with arbitrary arg(y4), and we take 
Im(/i) = 0. In part of SUSY parameter space we find enhancement of these cross sections, 
and hence an increase in direct and indirect detection rates; while in other parts of parameter 
space the cross sections are suppressed. 

The phase of A enters into the neutralino cross sections in two places: 1) into the squark 
masses, and 2) into the Higgs sector. For example, one of the processes that contributes 
to neutralino annihilation is s-channel exchange of the three neutral Higgs bosons, h, H, 
and A, into final state fermions. (see fig. 1). The first two of these neutral Higgs bosons, h 
and H, are CP even, while A is CP odd. The new aspect considered here is the possibility 
of mixing between the CP-even Higgs scalars {h and H) and the CP-odd scalar A. This 
mixing was first studied by Pilaftsis who found that the size of CP violation can be 
fairly large, i.e. of order one, for a range of kinematic parameters preferred by SUSY. He 
found that a large HA mixing can naturally occur within two-Higgs doublet models either 
at the tree level, if one adds softly CP-violating breaking terms to the Higgs potential, or 
at one loop, after integrating out heavy degrees of freedom that break the CP invariance of 
the Higgs sector, such as heavy Majorana neutrinos. In any case, in this paper, we consider 
the one-loop effects of Im(y4) 7^ on scattering and annihilation cross sections relevant to 
direct and indirect detection. 

In Section II we discuss our general approach. In Section HI, we discuss the squark sector, 
and in Section IV, the Higgs sector. In Section V, we discuss experimental constraints on 
the parameters. We present our results in Section VI. 

II. GENERAL APPROACH 

The minimal supersymmetric standard model provides a well-defined calculational frame- 
work P, but contains at least 106 yet-unmeasured parameters 0. Most of them control 
details of the squark and slepton sectors, and can safely be disregarded in dark matter 
studies. So similarly to Bergstrom and Gondolo we restrict the number of parameters 
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to 6 plus one CP violating phase: the "CP-odd" scalar mass tua (which in our CP vio- 
lating scenario is just a mass parameter), the Higgs mass parameter fj,, the gaugino mass 
parameter M2 (we impose gaugino mass unification), the ratio of Higgs vacuum expectation 
values tan /?, a sfermion mass parameter M (not to be confused with the sfermion mass, see 
eqs. (0) and (^ below), and a complex sfermion mixing parameter A = At = At, for the 
third generation (we set the A's of the first two generations to zero). The phase of A is the 
only CP violating phase we introduce besides the standard model CKM phase. 

We use the database of points in parameter space built in refs. [pHTOD, setting their A^ 
equal to Af. Hence we explore a substantial fraction of the supersymmetry parameter space, 
running through different possible neutralinos as the lightest SUSY particle. 

We modify the squark and Higgs couplings in the neutralino dark matter code DarkSUSY 



T^ to include a non-zero phase of A. We also add all diagrams that contribute to neutralino 



scattering and annihilation and would vanish when CP is conserved. 

To investigate the effects of the phase of A, we perform the following procedure. For 
each of the 132,887 sets of parameter values in the database, we run through 50 values of 
the phase of A, so that we effectively explore 50 x 132, 887 ~ 6.6 x 10^ models. We loop over 
a circle with arg(^) varying from to 271. At each point, we check bounds on the electric 
dipole moment, on the Higgs mass, on other particle masses, on the 6 — >■ 57 branching ratio, 
and on the invisible Z width (table I gives a listing of the bounds we apply). If any of these 
bounds are violated, we move to the next point on the circle. If all the bounds are satisfied, 
we calculate the spin-independent neutralino-proton scattering cross section a^p. We record 
the two values of the phase of A where a^p is highest and lowest, respectively, with the 
bounds satisfied. Then, once we have looped through all the possible values for the phase 
of A, we have found the two points with the maximum enhancement and suppression of the 
scattering cross section. We then compare with the scattering cross section in the case of 
no CP violation. 

We do the same for the annihilation cross section times relative velocity av at relative 
velocity v = (we recall that (T~l/fasf— s-0). Thus we obtain the values of the phase of 
A where av is maximum and minimum. 



III. SQUARK SECTOR 

The (complex) scalar top and bottom mass matrices can be expressed in the {qhAR) 
basis as 

/^| + '^' + (^3g-e,sin2^H/)cos2/3m| (a; - /ii?,) \ 

^ V (^g - l^*Rq) M^ + ml + eg sin^ Owml cos 2/3 / ' 

where g = t or 6; ej = |; Cf, = — |; T-^t = \] T^b = Rt = cot (3, Rb = tan/5; and M~ = M~ 
[M|] for t [b]. We set 

= = M5 = M, (3) 

our sfermion mass parameter. Even in the case of no CP violation, when both /i and A are 
real, there is mixing between the squarks, and this matrix must be diagonalized to find the 
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mass eigenstates. Here we take A to be complex. Then we obtain the mass eigenstates gi, q2 
from the weak eigenstates q^, through the rotation 



cos 6g sin OgC^^'i 
— sin 9g cos 6'ge*'^« 



(4) 



where 7^ = arg(y4* — fiRq) and the rotation angle 9q (— 7r/4 < 6q < may be obtained 
by 



tan(2^q 



2'mq\A\- [lRq\ 



MI-MI + {2eq sin^ Ow - T^qWz cos 2/? ' 



(5) 



The masses of qi and q2 are then given by 
ml „ = ^ i M| + M| + T-iqml cos 2(3 



11,2 



± sign(0,-) ^/ M| - + (2e, sin^ Ow - T^qWz cos 2/3 + Ami - /^^' 



(6) 



The + sign is for q\ and the — sign for g2- 

The mixing in eq. (H) also modifies the squark couplings to the neutralino and the 
corresponding quark. Writing the relevant interaction term as 

(7) 



(8) 



>Cint = qiX [gixQ^L + dlkxQ^R) Q + h.c, 
with Pl = {1- 75)/2, Pr = (1 + 75)/2, and i = l,2, we have 



K 



COS Qq sin Qqt^^^ 
— sin Qq cos 6'q;e*'''« 



where K = L,R, 

9LL = -V2 {TsqgN,2 + {cq - T3q)g'Nu) , 

and 

gm^Nu 



9kl\ 

gKRj 



gRR = V2eqg'Nn, 



gLR — gRL 



\^m\Y sin P 



(9) 



(10) 



for the up-type quarks, 



gLR — gRL 



grUdN: 



13 



\/2mw cos P 



(11) 



for the down-type quarks. 

The expressions in this section apply to sleptons provided up- type (s) quarks is replaced 
with (s)neutrinos and down-type (s)quarks with charged (s)leptons. 
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IV. HIGGS SECTOR 



A. Higgs masses 

We evaluate the Higgs boson masses in the effective potential approach. The radiatively 
corrected Higgs boson mass matrix can be written as 



f m% cos^ 13 + m\ sin^ /3 + An — {m\ + m|) sin /3 cos /? + A12 A13 
-(m^ + m|) sin/?cos/? + A21 m| sin^ /5 + cos^ /? + A22 A 
A31 A32 



23 



:i2) 



in the basis Hi, H2, H^. Here Ajj = Ajj are the radiative corrections coming from quark 
and squark loops, with A13 and A23 arising from CP violation. We take An, A12, and A22 
from ref . . 



An 



^22 



^12 



^9' 



4 

cos^ [3 
-\ 



In- 



2 2 



b2 



+ 2Zh\n 



2 

mi 

01 



ml 



m 



b2, 

4 



m 



cos^ P 
4 / 

In 



sin P 



m 



2 2 
mf mj 

mf 



2 ' 
ti 



sin^/? sin=^/5 ''1' 



sin P 



t Wt iln"^ + Z,g{ml,ml: 



m 



t2 



ml 



m 



W, \ In -^ + Z,g{mlmi 



cos^ P 



m 



62 



(13) 



(14) 



(15) 



where 



g{ml,ml 



Re{fiAq) — |yUp_Rq 
2 2 ' 

\Ag\^ -Re{i^Ag)Rq 



m 



92 



91 



2 + m2 J ml 

mf — m2 ml 



(16) 
(17) 
(18) 



We have rewritten A13 and A23 from ref. |^ in a way that shows their proportionality to 
lm{iiA). 



A 



fe3 



gAqiqi { 77 {gHkqLQL + fi'^^i 



k^RlR 



log 



91 



'72 



+ 



sin2^<jRe(e'>^/^,g^gJ + -cos2eg (^/f.^^gi - gn.qam) 



(19) 
(20) 
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where the couphngs of the Higgs bosons to the squarks are 

^ grnj Im(//^4f) 

9Ahh ™ „• 2 o ^2 ™2 ' 

gml Im(/iAfe) 
mw cos^ B m~ — m~ ' 

01 02 

= 7— I Tst - sm 6'^/ cos /3, 



dHltLtL 
tRtR 



cos 6^1^ 

= — et sm Ow cos p, 

cos 

grutfJ.* 



n T- 1^36 - Cfe sm Ow ) cos /3, 

rriTfyr cos p cos o'vy ^ ^ 

,-,2 



SHlbLbL ~ 



gml gmz . a 
9HibRbR^ 75 sm 9w cos 

^ gm^ gmz 



(Tst - et sin^ 6w) sin /3, 



9H2tRtR — 



gml 



mw cos p cos t^vt^ 



_ gmtAt 
gH^tRh- 2mtysin/3' 

9H2brb, = Z~ i^3b - Cfe sm Ow j sm /3, 

COS^H/ ^ ^ 



(21) 

(22) 

(23) 
(24) 

(25) 

(26) 

(27) 

(28) 
(29) 

(30) 

(31) 
(32) 
(33) 
(34) 



9hJ^rH 2mwC0sP' ' ' 

Neglecting D terms, as we should for consistency with the CP even part and the vertices 
in our effective potential approach, the corrections A13 and A23 simplify to 



Ai3 = 



^g' 



cos^ P 



■6,. "fa 



sin^/3 



(35) 



A23 = 



sin=^/3 



??7. 



(36) 



with 
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^, = (37) 



ml - 



The key thing to notice is that the self-energies are proportional to lm{nA). For n 
real, they are hence proportional to Im(y4). 

We use the effective potential approach to obtain the Higgs masses and couplings. The 
Higgs mass eigenstates hi {i = 1, 2, 3) are obtained by diagonalizing the Higgs mass matrix 
including radiative corrections in eq. (^) through the orthogonal Higgs mixing matrix O as 

Hi = Oijhj (38) 

In practice, it is convenient to implement the diagonalization in two steps, to separate the 
CP violating contributions. First we diagonalize the "CP-even" part through 

= 0° (39) 

where = H,h,A for z = 1, 2, 3 respectively. The matrix would be the Higgs mixing 
matrix in absence of CP violation 



0° = I sina cosa | , (40) 
with 

2 A/f2 

Then we further rotate to the mass eigenstates with an orthogonal matrix O' as 

= Ol^h, (42) 

with O' = 00^^. This two step procedure allows for a rapid introduction of CP violating 
mixing angles for the Higgs sector in the DarkSUSY code. 




B. Higgs couplings 

We will include CP violating effects by rotating couplings of Higgs particles to other 
particles as described in this section. In the effective potential approach we neglect ver- 
tex corrections. This incorporates the dominant corrections of 0{g'^m^ /m^r), and neglects 
corrections of C((y'^m^/m^). 

There are terms in the Lagrangian that couple the Higgs particles to other particles that 
are linear in the Higgs fields, for example 

g^iqq^m = 9h,qqOjihjqq- (43) 

Terms of this type include coupling to fermions, as shown above, and also terms such as 
gwH+^iWH~^^i. We will define rotated couplings via 
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9h,ab = (44) 

where a and b stand for the appropriate particle name. 
Those terms with two Higgs bosons in them, such as 

gz<s>3^,Z^3d^^i = gz<s,,<s,p'k30jiZhkd^hj, (45) 

must have the couphngs rotated with two multiphcations by O', e.g., 

gzHh, = gzA'i>P'jiO'k3 -{k^ j). (46) 

Note that, in this particular term, the appropriate antisymmetry properties are maintained, 
and i takes on values 1 or 2 only. 

We have carefully rotated all couplings involving one, two, or three Higgs bosons. It 
is these rotated couplings that we use in the numerical code. (i.e. we replace the ordinary 
Higgs couplings with these rotated couplings.) 

As an example, we give the Higgs-quark and Higgs-neutralino vertices that appear in 
the neutralino-proton spin-independent cross section. 

ghiuu = (0-1 sin a + cos a + 0-3^ cos/?), (47) 

ZTn\Y sm p 

gh^dd = - 2m^cos(3 " ~ *^»2 sin a + O-^i sin /3) , (48) 

gh^xmxn = ^igN^2 - g'N*mi) [K^ {-0[^ cosa + 0-2 sin a + O'J sin/?) + 

A^*4 {p[i sin a + 0^2 cosa — O'^^i cos /?)] + (m ^ n). (49) 

Here u stands for down-type quarks and neutrinos, d stands for up-type quarks and charged 
leptons. 



V. EXPERIMENTAL BOUNDS 
A. Bounds on masses 

We impose experimental bounds on the invisible width of the boson, F^^, and on 
particle masses as listed in table I. 

Since the h, H, and A are rotated into new mass eigenstates bosons, we use the most 
model independent constraint on the neutral Higgs masses: we take rrih^ > 82.5 GeV. This 



constraint was reported by the ALEPH group [T^ at the 95% C.L. as a bound on all Higgs 
masses, independent of sin^ (/? — «). Note that this bound, which is a 10% improvement over 
previous bounds, renders the cross section for direct detection of SUSY particles smaller by 
a factor of two. This suppression arises because the dominant contribution to the scattering 
cross section is via Higgs exchange and scales as a^p oc 1/iti\_. 
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B. Bounds on CP violation 



57), and on the electric dipole 
with inclusion of the one-loop QCD 



We impose bounds on the branching ratio BR(6 
moments of the electron and of the neutron and 

For BR(6 S7) we use the expressions in ref. |]14 , 
corrections. 

Since we assume that the only new CP violating phase is that of A, the leading contri- 
bution to the electric dipole moment (EDM) arises at two-loops |]T5[. Chang, Keung, and 



Pilaftsis [^] have calculated two-loop contributions to the electric dipole moment (EDM) 
which originate from the potential CP violation due to a nonzero phase of A. We rewrite 
them showing explicitly their dependence on Im(//A). They find the electric and chromo- 
electric EDM of a light fermion / at the electroweak scale as 



(df. 



EW 



q=t,h 



(50) 



f JEW 



m. 



/ 2 



12 



(51) 



where aem = e^/(47r) is the electromagnetic fine structure constant, = gl/i^in) is the 
strong coupling constant, all the kinematic parameters must be evaluated at the electroweak 
scale m^, is the electric charge of particle i, Rf = tan/3 for / = u,c,t, Rf = cot /3 for 
f = e, fj,, T, d, s, b, and F{z) is a two-loop function given by 



Fiz) 



1 , x(l — x) 
dx- 



z — xil — x) 



In 



xil — x) 



(52) 



The EDM of the neutron can then be estimated by a naive dimensional analysis |[T6|JT7[] as 



dn = A{K-d^)+^<^^[Ad^,-d<^^). 



127r 

We take the numerical values r]^ = 1.53 and 77*" = 3.4 [p!7| . 
The CP violating quantities and are given by 



(53) 



6 



gm^ lm{fiAt 



and 



2m^ sin73(m| — m|) 

gml lm{fxAb) 
2m^ cos^/5(m?^ — m?J 



(54) 



(55) 



b2' 



As an upper bound to the contribution to the measured value of the electron EDM we take 



141 < 0.4 X lO-^^'ecm |§. The bound on the neutron EDM is \dn\ < 1.79 x lO'^^ecm p 
We keep only models that satisfy these bounds. 
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VI. SCATTERING CROSS SECTION 



The neutralino-proton scattering cross section for spin-independent interactions can be 
written as 



a 



XP 



TT 



(56) 



where n^^p = m^mp/ {m-^ + nip) is the reduced neutrahno-proton mass, and 

fqrnp 



G 



XP 



E 



q mq 



.1=1 hi 



1 ^^^9qkxq9^kxq^ 

^ fe=l ^^Qfc 



(57) 



We take f35 



The sum over q runs over all quarks. The coupling constants are given in eqs. (|8|,^7|-^9|) . 

/„ = 0.023, /d = 0.034, /. = 0.14, = = = 0.0595, (58) 



and 



m„ = 5.6MeV, = 9.9MeV, = 199MeV, 
rric = 1.35GeV, = 5GeV, rrit = 175GeV. 



(59) 



Notice that only the real part of the couplings of the Higgs and neutralinos to Higgs 
bosons in eq. (|^-^) enter the scattering cross section. Since both gAqq and Qaxx purely 
imaginary (because Im(/i) = 0), introducing a phase in A cannot possibly enhance the 
Higgs couplings in eq. (p7D. Similarly, the neutralino-squark-quark couplings can only be 
suppressed for Im(A) ^ 0. However, enhancements to the scattering cross section can still 
come from the Higgs or squark masses in the denominator in eq. (|5^). 



VII. ANNIHILATION CROSS SECTION 

The neutralino-neutralino annihilation cross section times relative velocity av is relevant 
for neutralino annihilations in the center of the Earth and Sun and in the galactic halo. An 
enhancement in av may lead to a higher annihilation signal from the Earth when the capture 
of neutralinos in the core has not yet reached equilibrium with their self-annihilation. An 
increased av gives directly an increased intensity of positron, antiproton, and gamma-ray 
fluxes from neutralino annihilation in the galactic halo. 

The neutralino annihilation cross section also determines the relic density of neutralinos. 
In this case, there are important contributions at w 7^ (p-waves, etc.) in large regions of 
the supersymmetric parameter space. Due to the excessive computational cost of obtaining 
the relic density in presence of CP violation, in this paper we consider only the v = case, 
and postpone the study of the effect of CP violating phases on the neutralino relic density. 
The enhancements and suppressions of at f = that we obtain in the following are 
indications of analogous enhancements and suppressions in the neutralino relic density. 

The annihilation cross section at f = includes the following contributions 
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E 



a- 



ff 



3 3 
i=l ij = i 



(60) 



where axY refers to the annihilation channel XX XY , which is open when 2m-^ > mx + 
my- 

The annihilation cross section in each channel can be written in terms of helicity ampli- 
tudes A as 

XxY 



axYV 



12871772^ helicities 



(61) 



where the amplitudes are normalized as in ref. pO| and 



A 



XY 



\ 



1 - 



{mx + mvY 



1 - 



{mx — myy 



4m^ 



(62) 



The DarkSUSY code already includes analytic expressions for each helicity amplitude re- 
quired in eq. (|6ll), with arbitrary complex couplings between the particles. Hence once we 
have rotated all vertices as described in sect. IV, and have added all annihilation diagrams 
that vanish when CP is conserved (e.g. the s-channel exchange of all Higgs bosons), the 
annihilation cross section including CP violation is automatically calculated correctly by 
DarkSUSY. 

For future reference, we list the individual contributions to the annihilation cross section 
including terms that violate CP. 



327r 



E 



4Im(^h,//)Im(c/;,^x,^J 4g^ffRe{gzxixi)i^f/^> 



^ ml^ -Aml-imhTh^ 



+ 



mi. 



+ 



+ E 



z 

« 12 + \glJ^){mf/m^) + 2Re{g'^^^^gf* 



fsXf' 



fsXf fsXf 



s=l 



mf + m 

Js 



+ 



NjXj^ml 



327r 
Xww 



j2 ^^'^^i9h,ff)'^'^i9h,xi 



X 



XI- 



m J — im F 



+ 



+ E 



^ - Am^^ - imhj^ m'j^ + m"^ - mj - im J 



+ 



(TZZV 



+ 



«7r 
Xww 

Xww 
167r 

Xww 



Xwwm^{\g^^^ 



+ \9- 



Wxxi 



P) + 2zm,.Im«^,,4;,,) 



c=l 



E 

c=l 
4 

E 



m_+ + mi 

Xc > 

[2{m^/mwf - \m{g 



m 



w 



im-+T-+ 

Ac Ac 



+ 



WxxI^'Wxx 



mi+ + m: 

Xc ^ 



m 



w 



2m-+F-4 

Ac Ac 



4 2Xzzm^\gzxxS - 2^"^xn Im(c/|^^„ 



m^^ -|- ml 



mr 



^mxrXxr. 



n=l ""Xn ' ""X 

4 ro/'^ /^...^2 11^ Tv^/'^2 



+ 



E 



[2{mjmwr - l]"^xn Im(c/|^^^ 



Stt 

(Th+w-v = (Jh-w+v 



n=l 



m?, + mi — m 

X" A 



^\±w^'x 



i^Xn'^Xn 



IGnm'^ 



w 



E 

i=l 



^igwKH± Im(5fh 



4m2 - imhTh, 



(63) 



(64) 



(65) 
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+ E 



(„R „R* 



^Wxxt^Hxxi 



+ {g- 



R* 

Wxxt^Hxxt 



SHxxt^H^xt^'^xt 



+ mi 



(m 



im-+r-+ 

A.C Ac 



E 

ra=l 



-2 Re{gzxx9*h. 



iXiXn 



]m^ + 2 Re{gzxx9h iXlXnj'^Xn 



+ mi 



K. +m|)/2 



^ 4:igzh,h, 1^9 h 



jXX)^X 



Zxx) 



A 



0-h,hjV 



hA 



GArcml 

X 



71=1 



- Ami - imij^h 



mi 



mr 



XlXn9hjXlXn)^X^Xr 



mt 



iXlXnghjXlXn, 



mt 



+ mi 



[ml + ml)/2 



Xn Xn 



+ E 



kXX)^X 



2 

, mi 

1 "fc 



4m?, 



+ 



igzh,h,{mj - m]) Re{gzxx) 



mr 



(66) 



(67) 



(68) 



Nf is 3 for quarks and 1 for leptons, ghiff and ghiXmXn given in eqs. (|47|-^9D, gj^^^ and 
gf are given in eqs. (p|-[Tl|), and 



'fxf 



gwKH± = 7T [O'ii sin(a - (3) + 0-2 cos(a - (3) + iO[^^ 



3J 5 



gZhihj 



9zff 



g 



Wxxt 



^Wxxi 



2 cos 9w 
2 cos 9\Y ' 



[O^]^ sin(a — /5) + 0-2 cos(a — 



V2 



+ gNi2v:,, 
+ gN^2Ucu 



9ZXr, 



g 



2 COS 9w 



(69) 
(70) 
(71) 
(72) 

(73) 
(74) 



Here V and U are the chargino mixing matrices. 

For real yU and real gaugino masses, as we assume here, the terms containing 
^^^^wxxt^wxxt^ Im(^|^^J in eqs. (|4|-||) vanish. 

Notice that in the annihilation into fermion pairs, in the first terms under absolute values 
in eq. (^) (see fig. 1(a)), there can be contributions from all Higgs bosons hi for which the 
imaginary part of ghixx non-zero. Examining eq. (|49|) for the couplings, recalling that the 
matrix elements are real and that for real /i and real gaugino masses NuNij are also 
real, we see that the hi contributes when O^g is non-zero. In the CP-conserving case, this 
happens only for 2 = 3, i.e. for the A boson, while with CP violation this occurs also for 
i = 1 and i = 2. The annihilation into // then proceeds through exchange of all Higgs 
bosons, raising the possibility of resonant annihilation when 2m^ is approximately equal to 
the mass of any Higgs boson. This phenomenon is peculiar to CP violation. An example is 
given in fig. 9 below. 



12 



VIII. RESULTS 



A. Results for the elastic scattering cross section 

In Figure 2 we show the neutrahno-proton elastic scattering cross section as a function 
of neutrahno mass for the ~ 10^ values in SUSY parameter space that we consider. There 
is no CP violation in the lower panel (ImA = 0), while CP violation is allowed in the upper 
panel. Also shown are the present experimental bounds from the DAMA [RB] and CDMS 



collaborations as well as the future reach of the CDMS (Soudan) |3^, CRESST |^ and 



GENIUS experiments. In the upper panel, it is the maximally enhanced cross section 
(as a function of arg(y4)) that is plotted. The red (dark) points refer to those values of 
parameter space which have the maximum value of the cross section for nonzero Im(y4) and 
which are experimentally excluded at zero Im(A). The blue (grey) region refers to those 
values of parameter space which are enhanced when CP violation is included and which 
are allowed also at zero lm{A). The green (light grey) empty squares refer to those values 
of parameter space which have no enhancement when CP violation is included. From the 
existence of the red points we conclude that there are indeed points in SUSY parameter 
space which are ruled out experimentally when CP is conserved but are allowed when CP 
is violated. 

By comparing corresponding points in the upper and lower panels of figure 2, we notice 
that there can be enhancement or suppression of the cross section when we allow for CP 
violation. There are two types of enhancement: one in which the model without CP violation 
is allowed and another in which it is experimentally ruled out. In the first case, it is possible 
to define a ratio between enhanced and unenhanced cross sections, -Rmax = (^^^/(^max- 
the second case, when both cr(0) and (t{tt) are excluded, it is not possible to define the 
previous ratio. Here cr™^^ is the maximally enhanced cross section as one goes through the 
phase of A, and ajij^^^ = max[(T(0), <j{ii)] is the larger of the unenhanced CP conserving cross 
sections. We plot i?™^^ in figure 3 as a function of o"[]^ax- those models in parameter space 
that we have considered, we notice that the enhancement due to CP violation is at most a 
factor of two. 

In figure 3, we have also plotted the ratio -Rmin = '^™™/'^min' which is a measure of the 
maximal suppression of the neutralino-proton cross section when CP violation is included. 
Here cr™™ is the maximally suppressed cross section as one goes through the phase of A, and 
crj^in = min[o"(0), (7(7?)] is the smaller of the unenhanced CP conserving cross sections. We 
see that significant suppression of the scattering cross section, as low as 10""^, is possible. 

In figure 4, we show the dependence of these enhancement and suppression factors -Rmax 
and -Rmin on the phase 0^ of A. The points are plotted at those values of 0a at which the 
maximum or minimum of the scattering cross section occurs. 



B. Results for the annihilation cross section 

We also have obtained values for the neutralino annihilation cross section av for the case 
of CP violation through the phase of A. In figure 5 we show the maximum value of av 
obtained as we vary (pA as a function of neutralino mass. As in the analogous Figure 2 for 
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the scattering cross section, the upper panel inchides CP violation while the lower one does 
not. The distinction between red (dark), blue (grey), and green (light grey) points is as in 
Figure 2. 

Figure 6 shows the enhancement of the annihilation cross section via the ratio -R^^ = 
(av)^^/ (crv)^^. We see that the annihilation cross section can be significantly enhanced 
for CP violation with Im(yl) 7^ 0, by as much as a factor of 10^. In all models for which 
we find an enhancement in the annihilation cross section of at least 10^, the enhancement 
is due to an s-channel resonance with the exchange of one of the Higgs bosons hi, /i2 or h^. 
See fig. 9 for an example. 

A similar ratio can be constructed for to show that the suppression due to CP vio- 
lation can be roughly a factor of 50. The dependence of these enhancement and suppression 
factors and i?min on (pA are plotted in figure 7. 



C. Phase dependence of the results 

In the four panels in each of the fig. 8-11 we display the behavior of the scattering cross 
section a^p, the annihilation cross section av, the branching ratio BR(6 — > S7), and the 
lightest Higgs boson mass m/jj as a function of the phase (pA of A. In the third and fourth 
panels we hatch the regions currently ruled out by accelerator experiments. In all four panels 
we denote the part of the curves that is experimentally allowed by thickened solid lines, and 
the part that is experimentally ruled out (as seen e.g. in the third and fourth panels) by 
thinner solid lines. 

For the models shown in figs. 8-11, we give the values of the input parameters and of 
the neutralino mass and composition (gaugino fraction jiVnp + |A'^i2p) in table 2. 

In the case plotted in fig. 8, the possible phases are bound by the limit on the 6 57 
branching ratio. In the allowed regions, the scattering cross section at CP-violating phases is 
suppressed, while the annihilation cross section is enhanced. The latter takes its maximum 
allowed value when the 6 — > 57 limit is reached. 

Figure 9 presents another case in which the phase of A is bounded by the b s'-f branch- 
ing ratio. Here the scattering cross section is enhanced by only 2%, while the annihilation 
cross section is enhanced by a factor of ~ 222 at (pA — 0.1297r. This is due to a resonant 
annihilation of the neutrahnos through s-channel exchange of the hi Higgs boson (fig. 1(a)), 
which occurs when 2m^ — rrih^ (see the lowest panel). Notice that in the CP conserving 
case, the s-channel exchange of the CP-even hi boson vanishes at f = because for real 
XXhi couplings the amplitude is proportional to XX which is zero at f = 0. In presence 
of CP violation, the xxhi couplings are in general complex, and the amplitude contains a 
contribution from xibX which does not vanish at v — 0. So the hi resonant annihilation 
seen in fig. 9 is only possible when CP is violated. 

Figure 10 shows a case which is experimentally allowed for all values of the phase (pA- 
The maximum of the scattering cross section takes place at the CP-conserving value (pA = tt 
and the minimum at 0^ = 0. The annihilation cross section on the other hand is enhanced 
by CP violation, as can be seen in the second panel. Notice that its maximum occurs at 
(pA — 0.4l7r, which is not the point of maximal CP violation 0^ = 7r/2. 
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Finally fig. 11 displays an example in which both CP conserving cases are experimentally 
excluded while some CP violating cases are allowed. This is one of the red (dark) points in 
fig. 5. The 0^ = case is ruled out by the bounds on both BR(6 — 57) and the Higgs mass, 
the 0^ = vr case by only the bound on the Higgs mass. Notice that the scattering cross 
section is of the order of 10~^ pb, in the region probed by the direct detection experiments. 
The annihilation cross section peaks at (pA = 37r/4; notice that again this value is not the 
point of maximal CP violation (pA = 7r/2. 

IX. CONCLUSIONS 

We have examined the effect of CP violation on the neutralino annihilation and scatter- 
ing cross sections, which are of importance in calculations of the neutralino relic density and 
of the predicted rates for direct and indirect searches of neutralino dark matter. Specifically 
we have considered the case in which the only CP violating phase in addition to the standard 
model CKM phase is in the complex soft trilinear scalar couplings A of the third genera- 
tion. This phase affects the squark masses and through radiative corrections generates a 
mixing between CP-even and CP-odd Higgs bosons. This mixing modifies the neutralino 
annihilation and scattering cross sections in the kinematic regimes relevant for dark matter 
detection. Exploring ~ 10^ points in supersymmetric parameter space with a non-zero phase 
of A, we have found that: (1) the scattering cross section is generally suppressed, even by 7 
orders of magnitude in special cases; (2) the annihilation cross section can be enhanced by 
factors of 10^ as resonant neutralino annihilation through a Higgs boson becomes possible 
at CP- violating values of the phase of A. We have also found cases which are experimen- 
tally excluded when CP conservation is imposed but are allowed when CP conservation is 
violated. Some of these cases have neutralino masses and cross sections in the region probed 
by current dark matter searches. 
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(a) (b) 

FIG. 1. Processes that contribute to neutralino annihilation and scattering and are affected 
by the CP violating phase of A. For annihilation: (a) s-channel diagrams via the three neutral 
Higgs bosons hi,h2, and hs into final state fermions /, and (b) t-channel diagrams via intermediate 
squarks / into final state fermions. For scattering: crossed diagrams. 
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Bound Ref. 
rf^ < 502.4MeV @" 



mH± > 59.5GeV pT 



mh, > 82.5GeV p| 



m-+ > 91GeV if m^o - m-+ > 4GeV pi 

m -+ > 64GeV if m^o > 43GeV and m -+ > m^o [|3| 

m-+ > 47GeV if m^o > 41GeV ' H 

m-+ > 99GeV |p5| 
2^ 



m^o > 23GeV if tan /3 > 3 H 
m^o > 20GeV if tan /3 > 2 @ 
m^o > 12.8GeV if mc < 200GeV H 



m^o > 10.9GeV 
m^o > 44GeV 
m^o > 102GeV 
m^o > 127GeV 

A4 



> 212GeV if mg, < |0| 

> 162GeV 111 



mg, > 90GeV if m-g < 410GeV ^ 

mg, > 176GeV if rUg < 300GeV H 

mg, > 224GeV if m-g > m-g H 

me > 78GeV if m^o < 73GeV |4| 

m^ > 71GeV if m^o < 66GeV [H 

mf > 65GeV if m^l < 55GeV || 

> 44.4GeV ' Pl 



1 X 10"^ < BR(6 ^ S7) < 4 X 10"^ 


m 


141 < 0.4 X IQ-^^ecm. 


11 


\dn\ < 1.79 X 10^25 gpjj^ 


11 



TABLE I. Experimental bounds we use in this paper. We do not include cosmological bounds 
nor bounds from dark matter searches. 



19 





JEsp4_001509 


JE27__004174 


JE28__002656 


JEsp4_002809 




Fig. 8 


Fig. 9 


Fig. 10 


Fig. 11 


^l [GeV] 


-331.433 


-271.973 


-234.128 


958.213 


M2 [GeV] 


390.064 


106.141 


338.688 


-153.256 


ruA [GeV] 


84.2527 


168.935 


325.691 


106.804 


tan/3 


31.6126 


4.37629 


1.80096 


48.4750 


M [GcV] 


1085.05 


494.379 


1856.43 


890.647 


A/M 


2.71920 


0.661158 


1.88819 


-2.05440 


[GeV] 


191.46 


54.95 


172.4 


77.04 


|iVii|2 + |iVi2|2 


0.9459 


0.9806 


0.9571 


0.99786 



TABLE II. Model parameters and neutralino mass and composition (gaugino fraction) for the 
examples in figs. 8-11. 
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10 102 IQ3 IQ4 



m 
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GeV 



FIG. 2. Neutralino elastic scattering cross section (in pb) as a function of neutralino mass (in 
GeV) for 10^ values in SUSY parameter space. The upper panel is for the case of CP violation 
via Im(A) ^ while the lower panel is for the case of no GP violation. In the upper panel, it is 
the maximally enhanced cross section (as a function of arg(A)) that is plotted. The red (dark) 
points refer to those values of parameter space which have the maximum value of the cross section 
for nonzero Im(A) and which are experimentally excluded at zero Im(A). The blue (grey) region 
refer to those values of parameter space which are enhanced when GP violation is included and 
which are allowed also at zero Im(^). The green (light grey) empty squares refer to those values 
of parameter space which have no enhancement when CP violation is included. The solid lines 
indicate the current experimental bounds placed by DAMA and CDMS; the dashed lines indicate 
the future reach of the CDMS (Soudan), GENIUS, and CRESST proposals. 
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FIG. 3. Enhancement and suppression of elastic scattering cross section for the case of CP 
violating arg(A). The plot shows the ratio i?max = cr^'^^/Cmax ^ function of the unenhanced 
scattering cross section a^^^^ = max[iT(0), cr(7r)]. Here fj™'^^ is the enhanced scattering cross section 
and the superscript max indicates the maximal enhancement as one goes through the phase of A. 
The denominator of the ratio -Rmax chooses the larger value of the scattering cross section without 
CP violation, i.e., for phase = or phase = vr. Similarly, the ratio i?min = '7™™/'^min plotted; 
this time the denominator chooses the smaller value of the scattering cross section without CP 
violation, (tJ^jj^ = min[cr(0), 0"(7r)]. 
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FIG. 4. The enhancement /suppression factors i?max and i?min defined in the caption of figure 
3 as a function of the values (j)A of the phase of A where the maximum/minimum occur. 
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FIG. 5. Same as fig. 2 but for the neutralino annihilatfon cross section times relative velocity av 
(in cm^/s at f = 0) as a function of neutralino mass (in GeV) for ~ 10® values in SUSY parameter 
space. 
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FIG. 6. Enhancement and suppression of neutralino annihilation cross section for the case of 
CP violating arg(^). The plot shows the ratio R^^^ = (f'"'y)™'^^/(c^^)max a function of the 
unenhanced annihilation cross section {crv)^^^^ = max[(Tu(0), cru(7r)]. Here (av)'^'^^ is the enhanced 
scattering cross section and the superscript max indicates the maximal enhancement as one goes 
through the phase of A. The denominator of the ratio -R^ax chooses the larger value of the 
scattering cross section without CP violation, i.e., for phase = or phase = vr. Similarly, the ratio 
R^^ = (o-t;)™''^/(o-i;)|^i^ is plotted. Here iav)'^^^^ = min[o-t;(0), C7?;(7r)]. 
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FIG. 7. The enhancement /suppression factors R^^^ ^^'^ ^min defined in the caption of figure 
6 as a function of the values (pA of the phase of A where the maximum/minimum occur. 
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FIG. 8. The four panels from top to bottom display the following: the scattering cross section 
a^p in pb, the annihilation cross section av in cm^/s, the branching ratio BR(6 57) x 10^, and 
the lightest Higgs boson mass nih^ in GeV as a function of the phase (pA of A. CP conserving 
phases are 0a = 0, tt while all other values are CP violating. In the third and fourth panels we 
hatch the regions currently ruled out by accelerator experiments. In all four panels we denote 
the part of the curves that is experimentally allowed by thickened solid lines, and the part that 
is experimentally ruled out by thinner solid lines. In this figure, the possible phases are bound 
by the limit on the b ^ sj branching ratio. In the allowed regions, the scattering cross section 
at CP-violating phases is suppressed, while the annihilation cross section is enhanced. The latter 
takes its maximum allowed value when the 6 — ^ 57 limit is reached. 
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FIG. 9. Same notation as fig. 8. The phase of A is bounded by the b ^ sj branching ratio, 
the scattering cross section is enhanced by only 2%, and the annihilation cross section is enhanced 
by a factor of ~ 222 at (/)a — 0.1297r, where the annihilation proceeds through the hi resonance at 
2m-^ = m/ij (see bottom panel). 
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FIG. 10. "The Duck. " Same notation as fig. 8. This case is experimentally allowed for all values 
of the phase ^a- The maximum of the scattering cross section takes place at the CP-conserving 
value (j)A = 71" and the minimum at (pA = 0. The annihilation cross section is enhanced by CP 
violation, as can be seen in the second panel. 
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FIG. 11. Same notation as fig. 8. Here, both CP conserving cases are experimentally excluded 
while some CP violating cases are allowed. This is one of the red (dark) points in fig. 5. The 
scattering cross section is of the order of 10~^ pb, and lies in the region being probed by direct 
detection experiments. The annihilation cross section peaks at (f)A = Sir/ 4; notice that this value 
is not the point of maximal CP violation (j)A = 7r/2. 
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